Real-time density matrix renormalization group dynamics of spin and charge 
transport in push-pull polyenes and related systems 
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In this paper we investigate the effect of terminal substituents on the dynamics of spin and charge 
transport in donor-acceptor substituted polyenes (D — {CII)x — A) chains, also known as push-pull 
polyenes. We employ a long-range correlated model Hamiltonian for the D — (CII)^ — A system, 
and time-dependent density matrix renormalization group technique for time propagating the wave 
packet obtained by injecting a hole at a terminal site, in the ground state of the system. Our 
studies reveal that the end groups do not affect spin and charge velocities in any significant way, 
but change the amount of charge transported. We have compared these push-pull systems with 
donor-acceptor substituted polymethine imine (PMI), D — {CHN)x — A, systems in which besides 
electron affinities, the nature of pz orbitals in conjugation also alternate from site to site. We note 
that spin and charge dynamics in the PMIs are very different from that observed in the case of 
push-pull polyenes, and within the time scale of our studies, transport of spin and charge leads to 
the formation of a "quasi-static" state. 

PACS numbers: 72.15.Nj, 72.80.Lc, 71.10.Fd 



INTRODUCTION 



Conjugated organic materials have a variety of ap- 
plications. They are already being used in organic 
light emitting diodes (OLEDS) and organic thin-film 
transistors pi"— They are also considered to be potential 
candidates for single-molecule electronic and spintronic 
devices. The issue of spin and charge transport in tt- 
conjugated organic systems, therefore, is of prime impor- 
tance. Thus, it is both important and necessary to under- 
stand the role of chemical modifications (substituents) on 
the dynamics of spin and charge transport in these sys- 
tems. Electrons in these molecular materials experience 
strong electron correlations by virtue of reduced dimen- 
sionality and these interactions are also long-ranged since 
the systems are semiconducting. Therefore, to study 
transport of charge and spin in these systems, it is nec- 
essary to have appropriate models as well as techniques. 
Dynamics in interacting systems can be studied either 
in the frequency or in the real-time domain. The dy- 
namics of strongly interacting systems in the frequency 
domain has been possible due to the Lanczos and correc- 
tion vector techniques^"— using both exact diagonaliza- 
tion method for small systems, and the density matrix 
renormalization group (DMRG) method for large sys- 
tems. 

Studying dynamics in the real-time domain throws a 
lot more light into the transport mechanism. Exact time 
dependent techniques are possible only for small inter- 
acting systems. The advent of the time-dependent den- 
sity matrix renormalization group (td-DMRG) method 
has allowed the study of out- of- equilibrium dynamics 
in large low-dimensional strongly correlated systemsi^"— 
Recently, we have investigated non-equilibrium dynamics 



of spin and charge transport in unsubstituted polyenes, 
which are a class of 7r-conjugated organic materials, us- 
ing the double time window targeting (DTWT) td-DMRG 
technique developed by usf^^ In the present paper we 
extend our studies to address non-equilibrium spin and 
charge dynamics in the push-pull polyenes, characterized 
by the presence of an electron donating (push) group, 
and an electron accepting (pull) group, schematically 
represented as D — {CH)x — A. Presence of the donor 
and acceptor groups polarizes the 7r-electron bridge such 
that there is partial charge transfer from the donor to 
the acceptor group. The electronic structure of a push- 
pull polyene can be described by two canonical valence 
bond (VB) structures, namely, a neutral polyene struc- 
ture and a zwitterionic structure, also known as charge 
transfer(CT) configuration, where an electron is trans- 
ferred from D to At^ This leads to the naive expecta- 
tion that these groups will have significant influence on 
spin-charge dynamics. The push-pull polyenes have been 
so far studied mostly in the context of nonlinear optical 
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In this paper we report our time-dependent DMRG 
studies on the spin and charge transport in push-pull 
polyenes and compare these results with those on poly- 
methine imines which are quasi one-dimensional systems 
with alternate C and N atoms in the conjugation back- 
bone. The organization of the paper is as follows: In the 
next section we provide details about the model Hamil- 
tonian and the computational method used in this study. 
In Sec. Ill we present our results with discussions. Sec- 
tion IV concludes the paper. 
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II. MODEL AND COMPUTATIONAL 
METHODOLOGY 



The Pariser-Parr-Pople (PPP) HamiltoniaiJ^ with 
dimerization and site energies is appropriate for describ- 
ing the low energy physics of 7r-conjugated systems. The 
PPP model Hamiltonian is given by 
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Here, N is the number of carbon atoms in the polyene 
chain, (^4,(7) creates (annihilates) an electron with 
spin orientation a in the Pz orbital of the i**^ carbon atom, 
to is the average transfer integral and, 5 {0 < 5 < 1) is 
the bond alternation parameter. The orbital energy and 
on-site Coulomb repulsion of the Pz orbital on the i*'^ 
carbon atom are given by and Ui, respectively and hi 
is the number operator on the i*'^ site. Vij is the intcr- 
site Coulomb repulsion between sites i and j, and Zi is 
the on-site chemical potential at the i*^ site. In case 
of unsubstituted polyene systems,— Ui = 11.26 eV, 
= and Zi = 1, for all sites, to = —2.4 eV and 5 = 
0.07. The intersitc interaction between electrons, Vij, is 
interpolated according to the Ohno schemCfii between U 
for ?■ = and — for r -> 00 as. 



14.397 



28.794 



-1/2 



(2) 



We have used single-bond length of 1.495 A, double- 
bond length of 1.299 A, and a bond angle of 120° be- 
tween successive bonds. These parameters have correctly 
predicted a whole range of properties of the low-lying 
states of conjugated systems in general and polyenes in 
particular.—"— When push and pull groups are intro- 
duced, we assume that only the energies of those Pz or- 
bitals to which the push and pull groups are attached, 
change due to inductive effect and all other parameters 
of the model remain unchanged. The donor group raises 
the orbital energy of the first carbon atom to which it 
is attached by -t-eu, while the acceptor group lowers the 
orbital energy of the last carbon atom to which it is at- 
tached by — ea, where ejj and e^i > 0. We assume that 
= —eA (symmetric push-pull polyenes). We have stud- 
ied push-pull polyenes of 30 and 40 carbon atoms with 
terminal push and pull groups, and have varied the push- 
pull strength |e|. Although, presence of the push and pull 
groups destroys both the electron-hole and the inversion 
symmetry, the total spin invariance of the Hamiltonian, 
remains preserved. 



For studying the dynamics of spin and charge trans- 
port, an initial wave packet | V'(O)) is constructed by 
annihilating an up spin electron from the first site of a 
push-pull polyene of N sites, in the ground state, | (/j^g). 



V^(O))-ci,t|0°). 



(3) 



The wave packet | "0(0)) evolves under the influence of the 
PPP Hamiltonian [Eq. (1)] following the time-dependent 
Schr"odinger equation and temporal dependence of site 
charge density {ni{t)), and site spin density (s|(t)), of 
this wave packet can be computed as. 



(4) 
(5) 



Here, | V'(^)) is the wave packet at time t. 

Real-time dynamics of the initial wave packet is stud- 
ied using the DTWT scheme)^ with the following specifi- 
cations: number of density matrix eigenvectors (DMEVs) 
retained, m ~ 300, time-step for evolution Ar ~ 0.0066 
fs, total evolution time T = 33.0 fs, the number of time 
steps in each window is kept at 130 which corresponds to 
0.858 fs, and the number of windows is 39. Although we 
compute charge (spin) density at all the sites, we focus 
only on the quantities (ni(t)) and (^^(t)) and, (si(t)) 
and (s|^(t)) at the terminal sites attached to the sub- 
stituents as these are sufficient for our purpose. 

Another class of donor-acceptor substituted systems 
examined are the polymethine imines (PMI), with the 
molecular formula D — {CHN)x — A.—"— These systems 
have alternately donor (C) and acceptor (N) atoms in 
conjugation. The bonding in this polymer corresponds 

to CH ^ N - CH ^ N and both carbon and 

nitrogen atoms are in sp^ hybridization. This system has 
been studied extensively for linear and non-linear optical 
properties. The nitrogen 2pz orbitals are lower in energy 
than the carbon 2pz orbitals and intra-orbital repulsion 
of the nitrogen 2pz orbitals are greater than that of the 
carbon 2pz orbitals, both these can be rationalized on 
the basis of the more compact 2p orbitals in nitrogen 
compared to carbon. The transfer integrals for C = TV 
and C - N are -2.767 eV and -2.317 eV, respectively; 
C = N and C-iV bond lengths are 1.273 Aand 1.425 A; 
Uc = 11.26 eV, Un = 12.34 eV; ec = 0.0 eV, = -2.96 
eV. The initial wave packet is constructed as before [see 
Eq. (3)] and time evolved by the DTWT scheme using a 
smaller time-step of Ar = 0.0033 fs, necessitated by the 
larger transfer integrals. 



III. RESULTS AND DISCUSSION 



In the absence of source and sink, the charge and spin 
of the hole decouple (spin-charge separation) and oscil- 
late between sites 1 and L, with time. Hence, the time 
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FIG. 1: (color online) Temporal variation of {ni{t)) and 
{si{t)) for dimerized {S = 0.07) push-pull polyenes. Panels 
(I) and (II) are for N = 30 chains while (III) and (IV) are 
for N — 40 chains. Panels (I) and (III) give evolution of 
charge while (II) and (IV) give spin density evolution. In 
each panel, bottom box gives the evolution of site 1 and the 
top box, evolution at the end of the chain L. In the charge 
density panels ((I) and (III)) |e| is 0.0 eV, 0.5 eV, 1.0 eV, 1.5 
eV and 2.0 eV for curves from top to bottom in the {ni{t)) 
boxes, and from bottom to top in the inL(t)) boxes. The 
same is the case in the spin density panels ((II) and (IV)). 
In color, |e| = 0.0 eV (black), 0.5 eV (red), 1.0 eV (green), 
1.5 eV (blue), 2.0 eV (maroon) curves. The positions of r^j^ 
and r£ in (I) and (III), and and r£ in (II) and (IV), 
are indicated with arrows. 



evolution profiles of (?^^(^)) f-nd (* = 1 and L), 

consist of a series of maxima and minima. The time t]J 
taken for the charge of the hole, and t£ for the spin of 
the hole, to propagate from site 1 to site L, is represented 
by the time at which the first major minimum (dip) ap- 
pears in the time evolution profile of {nL{t)) and (s|^(t)), 
respectively. Time required by the charge degree of free- 
dom, T^i^ , and spin degree of freedom, , of the hole to 
travel the round trip from site 1 ^ L — >■ 1 , is represented 
by the time at which the second major minima appears 



TABLE I: Variation in ^ and in the dimerized {5 = 0.07) 
Hubbard model with chain lengths 30 and 40, for ^ = 2.0, 
4.0 and 6.0. The same quantities for the 5 = 0.07 PPP model 
is also quoted. M and Q stand for model and quantities, 
respectively. 




^ = 2.0 ^ = 4.0 ^ = 6.0 PPP 



N 30 40 30 40 30 40 30 40 
i9i/i9F 1.33 1.32 1.35 1.39 1.30 1.18 5.05 4.80 
■&l/-dF 0.97 0.91 0.69 0.80 0.68 0.80 2.49 2.63 



in the time evolution profiles of {ni{t)) and (sf (t)). From 
these times, the charge and spin velocities of the injected 
hole are calculated as, d^l'^' — -4jh and -d^^j^' = -^fj^. 

The dynamics of spin and charge transport in systems 
with reduced dimensionality and electron-electron inter- 
actions, are expected to be different from each other, 
owing to spin-charge separation. When electron corre- 
lations are absent, the charge and spin velocities, and 
i?!^, are equal to the Fermi velocity, -dp. From our earlier 
td-DMRG studies^^ we have found that the Fermi veloc- 
ities of tight-binding chains of lengths 30 and 40 are 2.63 
and 2.72 (A/fs) respectively, for 5 = 0.0, and 2.43 and 
2.49 (A/fs), respectively, for 5 = 0.07. Table I presents a 
comparison of the ratio of ■d'l/^F and ^If-dp, for Hub- 
bard chains of 30 and 40 sites with different values 
and S = 0.07. These quantities are also computed for 
dimerized PPP chains with 30 and 40 sites, with |e| = 
0.0, and compared with those of Hubbard chains. As ev- 
ident from Table I, the values of the ratios of ■O'l/'dp and 
&1/'dp are much higher in the PPP model compared to 
the Hubbard model, due to the presence of long-range 
electron-electron interactions. 

As with unsubstituted polyenes in push-pull 
polyenes (Fig. 1) the spin propagates slower than the 
charge with t£/2l > t^/2l i'^l/2L < K/2l)- However, 
we find that the D and A groups have no significant effect 
on spin and charge velocities in the push-pull polyenes. 
Consequently, the charge and spin velocities in push-pull 
polyenes remain almost equal to those in unsubstituted 
polyenes (Table II). The push-pull substitution however, 
has an effect on the amount of charge transported, as 
seen from the depths of the minima in the time evolution 
profiles of {ni{t)) and {nL{t)). However, depth of the 
minima in (sf (i)) or {sf^{t)) do not change (see Fig. 1), 
implying that the push-pull groups have no influence on 
the spin degree of freedom of the hole. 

This implies that even though the spin and charge ve- 
locities in D — {CH)x — A polyenes remain largely un- 
affected by the donor-acceptor strengths, the amount of 
charge transferred in unit time decreases with increase 



TABLE II: Variation of t^; and r£ (fs) and, ■d'l and i9i (A/fs) 
with chain length (L), in symmetric push-pull polyene chains 
of length, 30 and 40 sites, with different values of |e| (eV). The 
PPP model parameters are, to — —2.4 eV, U = 11.26 eV, and 
5 = 0.07; L = 1.397(iV - 1 - S) A, N being the number of 
sites, and O stands for Observables. 
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40 
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Tl 



2.45 
4.96 
16.50 
8.15 



2.45 
4.92 
16.50 
8.22 



2.44 
4.90 
16.57 
8.25 



2.44 
4.87 
16.57 
8.30 



2.02 2.01 2.01 2.00 



3.35 
6.10 
16.24 
8.92 



3.34 
6.08 
16.28 
8.95 



3.30 
6.05 
16.42 
8.99 



3.25 
6.03 
16.74 
9.02 



2.42 
4.83 
16.71 
8.37 

2.00 

3.19 
6.01 
17.05 
9.05 



1.82 1.82 1.83 1.86 



in |e| as seen from the decrease in depth of the minima, 
when a hole is doped at the donor site. We can expect the 
opposite of when we dope an electron at the donor site; 
more charge will be transported from the donor in these 
cases as |e| increases. Our results indicate that push-pull 
substitutions do not alter the many-body character of the 
D — {CH)x — A polyenes, and that spin and charge veloc- 
ities are controlled by the correlation strength while the 
magnitude of charge transfer is controlled by the strength 
of the substitucnts. 

To understand in detail, the reason for the inability of 
the donor and acceptor groups to alter the many-body 
character of the D — {CH)x — A polyenes, we compute 
the inverse of the charge gap (AEc) and spin gap (AEs) 
of the push-pull polyenes. If Egs (A^ + 1) , Egg (A^ - 1) and 
Egs{N) denote the ground state energies of the {N + 1)- 
particle, {N — l)-particlc and A^-particle systems respec- 
tively, then the charge gap is given by, 

AE, = Eg,{N + 1) + Eg,{N - 1) - 2Eg,{N). (6) 

Similarly, the spin gap is defined as. 



AE, 



EoiSU = 3/2) - Eo{S^,, = 1/2), (7) 



where, Eo{S^^^ = 3/2) and Eo{S^^^ = 1/2) are the lowest 
energy states of the D — (CH)^ — A polyenes with S^^^^ 
3/2 and 1/2, respectively. In Fig. 2, we have plotted 

the variation in r^^' and 1/AE function of the 

strength of the push-pull groups. It is clearly observed 
that the donor-acceptor groups fail to alter either of the 
gaps, as a result of which, the velocities of the spin and 
charge remain unaffected. 

Comparing Figs. 1 and 3 it is clear that the dynam- 
ics of spin and charge transport in the D — {CHN)x — 
A chains is very different from that in the push-pull 
polyenes, owing to the different U values and site ener- 
gies of carbon and nitrogen atoms. For example, unlike 




FIG. 2: Variation in r^^" (fs), left-hand top and bottom plots, 
and l/Ai?c/s (eV~^), right-hand top and bottom plots, as a 
function of strength of push-pull groups (|e| = 0.0, 0.5, 1.0, 
1.5 and 2.0), for polyene chains of 30 and 40 sites. 



in push-pull polyenes (where the charge and spin of the 
injected hole after decoupling, oscillate back and forth be- 
tween sites 1 and L) in the D~{CHN)x~ A systems, the 
oscillations are damped and the spin and charge of the 
injected hole do not return after reaching the chain end. 
This is demonstrated by the observation that, although 
the first dip in the time evolution profiles of {nL{t)) and 
(s|^(i)) are well defined, there are no well defined sec- 
ond (major) minima in the {ni(t)) and (sf (t)) curves as 
a function of time. This may be because the hopping 
of an electron between unlike atoms is not a degener- 
ate resonant tunneling process and hence, is not strictly 
reversible. 

In order to investigate the fate of the injected hole, 
we study the temporal variation of total charge and spin 
densities in the left and right halves of PMI chains of 30 
and 40 sites, for |e| — 0.0 (Fig. 4). The total charge 
and spin densities for the left (right) half are defined as, 
{NLiR){t)) = EjeLiR) W), and (5£(^)(0) = E,eL(fl) 
{s'jit)) {L e [l,A^/2], R e [N/2 + l,iV]), respectively. 
These quantities serve as simple but effective probes to 
understand the motion of charge and spin of the hole, 
as they travel from the half of the system where charge 
injection occured to the opposite half of the system. If 
the time evolution profiles of these observables show sig- 
nificant oscillatory behavior, it implies that the charge 
and spin propagate back and forth between the ends of 
the system. However, if the oscillatory behavior is not 
pronounced and both the observables attain some aver- 
age value, it signifies that the charge (spin) moves in a 
such a manner that a "quasi-static" state is generated in 
which charge (spin) distribution of the system remains 
unchanged within the time of study. It is observed from 
Fig. 4 that for both the chain lengths, {S^{t)) goes from 
0.0 to —0.25 while, {Sf^{t)) goes from an initial value of 
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FIG. 3: (color online) Temporal variation in {ni{t)) and 
{si{t)) for polymethine imine system. Panels (I) and (II) 
are for = 30 chains while (III) and (IV) are for = 40 
chains. Panels (I) and (III) give evolution of charge while 
(II) and (IV) give spin density evolution. In each panel, 
bottom box gives the evolution of site 1 and the top box, evo- 
lution at the end of the chain L. In the charge density panels 
((I) and (III)) |e| is 0.0 eV, 0.5 eV, 1.0 eV, 1.5 eV and 2.0 eV 
for curves from top to bottom in the {ni{t)) boxes, and from 
bottom to top in the {nL{t)) boxes. The same is the case in 
the spin density panels ((II) and (IV)). In color, |e| — 0.0 
eV (black), 0.5 eV (red), 1.0 eV (green), 1.5 eV (blue), 2.0 eV 
(maroon) curves. In top boxes of the panels (I), (II), (III) 
and (IV), r£ and r[ have been indicated by black arrow. 



—0.5 to —0.25. Once these quantities have attamed the 
value of —0.25, they start to oscillate with time about 
this average value. The charge densities for the left and 
right halves of the PMI chains arc however, found to be- 
have little differently. In the case of the 30 site chain it 
is observed that, {NL{t)) and {Nii{t)) take the average 
values of ^^14. 25 and ~14.75, while for the longer chain of 
40 sites, these quantities oscillate about a mean value of 
'--'19.5. These data suggest that in the {CHN)^ systems, 
with increase in system size, within the computational 
time a "quasi-static" state is achieved more rapidly. 



FIG. 4: Variation in (A'^l{h)(^)) (plots in left box marked 
as (I)) and {S^(^jf^{i)) (plots in right box marked as (II)) 
with time, for unsubstituted PMI chains (|e| = 0.0) of 30 
and 40 sites. Solid curves correspond to {NL{t)} and {5'£(t)), 
respectively, and dashed curves correspond to {Nji{t)) and 



To understand this "quasi-static" state in more detail, 
we plot the time evolution profiles of charge and spin 
densities on carbon and nitrogen atoms in both halves 
of unsubstituted PMI chains of 30 and 40 sites (Fig. 5). 
These quantities are defined as, 



J6C 



JSC 



jew 



JEN 



(12) 

where L e [I, N/2] and R e [N/2 + 1, A^]. It is observed 
from Figs. 5(11) and 5(IV) that the oscillations in time, 
of {S'l / (j{t)) and (S^/j^ j^{t)), are opposite in phase, 
signifying that spin densities on the carbon and nitrogen 
atoms (in both half) are antifcrromagnetically coupled 
to each other. However, the total spin density carried 
by carbon atoms in both halves is found to be less than 
that carried by the nitrogen atoms. Furthermore it is 
seen that with time, both {Sf^ c(^)) ('^l Ar(0) de- 
crease, while (S'^ p(i)) and {S^ ff{t)) increase in mag- 
nitude, keeping the total spin density in either half at 
—0.25. From Figs. 5(1) and 5(111) we see that the av- 
erage value of {Ni^ip; (j{t)) is less than {Ni^/fij^{t)), and 
the total charge density carried by nitrogen atoms on the 
right half of PMI chains is more compared to that in the 
left half in line with the fact that the nitrogens are the 
"acceptor" atoms, being more electronegative. The op- 
posite is seen for the carbon atoms. The time evolution 
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time (fs) time (fs) 



(C/jv = 12.34 eV), the average site charge on N atoms 
is shghtly larger than 1.0 which is observed from the 
charge density distribution in the neutral ground state 
and the initial state. On the other hand, due to Uc = 
11.26 eV and ec ~ 0.0, the carbon atoms are slightly 
electron deficient in the {CHN)x system. Thus, propa- 
gation of the charge and spin degrees of freedom of the 
hole in this polarized background results in formation 
of the "quasi-static" state, which once formed, prevents 
both the charge and spin from returning to the site of in- 
jection as it requires reversal of polarization and is hence, 
energetically unfavorable. 

In unsubstituted PMl chains, t]J and t£ are 1.75 fs 
and 3.65 fs for N = 30, and 2.01 fs and 6.60 fs for N 
= 40. We find that the ■dh/'ds ratios for the TV = 30 
is 2.08 and for N = 40, 3.29. The 'dh/'&s ratio for a 
PMI chain of 30 sites agrees well with that of push-pull 
polyene with iV = 30 (|e| = 0). However, this ratio for 
= 40 (|e| = 0) is much higher in the PMI chain compared 
to push-pull polyene of same size. This is opposite to 
the behavior exhibited by push-pull polyenes, where the 
•dhj-ds ratio decreases with increasing chain length. In 
both polyenes and polymethine imines, the •dh/'&s ratios 
do not saturate, indicating that the 7r-coherence lengths 
for transport arc rather long. It also appears that al- 
ternate donor-acceptor sites along the chain enhance the 
velocities of both charge and spin. We have also studied 
substituents at the terminal sites in the {CHN)x systems 
and as with the polyenes, the strength of the push-pull 
groups do not affect the charge and spin velocities, or the 
nature of the "quasi-static" state. 



FIG. 5: Variation in total charge and spin densities of carbon 
and nitrogen atoms, with time, in the left and right halves 
of unsubstituted PMI chains (|e| = 0.0) with 30 and 40 sites. 
Top and bottom boxes correspond to A'' = 30 and 40 sites, 
respectively. Solid and dashed curves pertain to carbon and 
nitrogen, respectively. In the top box, plots marked (I) refer 
to charge density variation and plots marked (II) refer to 
spin density variation, in left (upper plots) and right (lower 
plots) half of the system, respectively. The same holds for 
plots marked as (III) and (IV), in the bottom box. 

profiles of total spin densities on carbon and nitrogen in 
both halves, compared to the total charge densities, are 
found to be more oscillatory. 

All these observations suggest that the "quasi-static" 
state is characterized by charge and spin of the hole dis- 
tributed between the left and right halves of the chain 
in almost equal amount, the positive or hole charge pre- 
dominantly residing on carbon atoms while the spin, on 
nitrogen atoms. The N atoms being more electronegative 
than the C atoms {cn = —2.96 eV, ec = 0.0 eV), prefer 
being electron rich by accommodating more charge den- 
sity. However, due to large on-site Coulomb repulsion 



IV. CONCLUSIONS 

To conclude, our studies show that the velocities of 
charge and spin transport are not affected by push-pull 
substituents, although the magnitude of charge trans- 
port depends on the push-pull strength. In the related 
{CHN)x system we find that, due to the underlying po- 
larized structure both the spin and charge transport are 
affected, and charge travels much faster than spin com- 
pared to push-pull systems. Furthermore, the transport 
of charge and spin results in the formation of a "quasi- 
static" state in which the injected hole resides in both 
halves of the system, almost in equal amount. We also 
note that in the PMIs the finite size effects are very large 
compared to the push-pull polyenes, suggesting much 
longer 7r-coherence lengths. 
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